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O 1 Abstract 

o. 

The determination of the spin distribution functions for the neutron requires an un- 
derstanding of the changes induced by a nuclear medium. We present the first estimates 
of the changes induced in the spin-dependent, valence parton distributions of the proton 
and neutron bound in 3 He and 6 Li. For 3 He the result is quite reassuring, in that the 
off-shell correction to g\ n is very small. On the other hand, in 6 Li the corrections to the 
spin distributions in the bound proton are sufficiently large that they should be taken into 
account if one aims for a precision of better than 10%. We provide a simple parametriza- 
tion of the off-shell changes in the valence spin distributions of 3 He and 6 Li, which can be 
used with any conventional estimate of nuclear binding and Fermi motion corrections. 



OO 



X 



There is considerable interest in deriving the corrections to the parton distributions of 
bound nucleons in order to finally understand the nuclear EMC effect [[l], [|. However, from 
the practical point of view there is a more pressing need. In the light of the activity associated 
with the violation of the Ellis- Jaffe sum rule 0, also observed first by the European Muon 
Collaboration (the "EMC spin effect") Q, there has been a concentrated drive to determine 
gi n in order to test the Bjorken sum rule ||. This has led to a tremendous increase in our 
knowledge of the proton and neutron spin structure, both theoretically and experimentally ||- 



TTfl . On the other hand, all of the information about the neutron comes from nuclear targets 
- so far D, 3 He and most recently 6 Li. It is therefore imperative to have control over the 
corrections to the quark distributions for a bound nucleon. 

For the deuteron the nuclear complications are sufficiently well understood that there 
have been several studies of the corrections to the neutron spin distributions extracted from 



D data [O. In this case it seems that the off-shell corrections are under control, at least for 
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x < 0.7. However, when it comes to heavier nuclei, such as 3 He, 6 Li and 15 N, the nuclear 
structure corrections may be of much more practical importance. Apart from this theoretical 
motivation, there is already some hint from recent experiments that nuclear corrections to 
nucleon properties may be significant in He. For example, recent studies of the neutron charge 
form factor in 3 He |l3j suggest a possible difference between the measurement in D and 3 He 
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In order to estimate the nuclear corrections to the parton distributions in nuclei such 
as these, at least at the present time when it is impossible to make a complete QCD based 
calculation, we need two ingredients. First, we need a model for the nucleon from which we can 
calculate free nucleon parton distributions that are in reasonable agreement with experimental 
data. Second, we need a theory based on the same quark model which can reproduce the 
essential features of the structure of the nucleus of interest. 

We consider first the issue of calculating the parton distributions corresponding to a 
particular quark model. There is a long history of studies of this kind, beginning with the 
work of Le Yaouanc et al. [I3|, Jaffe [|16| and Parisi and Petronzio |I7j in the mid-70's. The 
basic idea is to connect the twist- 2 parton distributions at a low scale with a valence dominated 
quark model - an idea exploited very successfully at the phenomenological level by Gliick et 
al. [13. A major problem with early quark model calculations, namely poor support, was 
solved about 10 years ago by formulating the problem in such a way that energy-momentum 
conservation was guaranteed before any approximation was made Using that technique 

it has been shown that, provided one allows for the hyperfine mass splitting between the S=0 
and S=l spectator pairs, the MIT bag model can give quite a good quantitative description of 



the observed valence parton distributions [20, 21 



The second essential ingredient has only been obtained quite recently. The quark meson 
coupling model (QMC), which is based on explicit quark degrees of freedom, is ideally suited 
for this purpose p2|-||27||. In its simplest form it is also based on the MIT bag model, with the 
interactions between different nucleons described by the exchange of scalar and vector mesons 
in mean field approximation. It has proven successful in reproducing the saturation properties 
of nuclear matter as well as the binding energies and charge densities of finite, closed shell 
nuclei. From the practical point of view, one of the most attractive features of the model is 
that it is not significantly more complicated than Quantum Hadrodynamics (QHD) |28| - even 
though the quark substructure of hadrons is explicitly implemented. A detailed description 
of the Lagrangian density, and the mean-field equations of motion needed to describe a finite 
nucleus, is given in Refs. [^3], p5fl . 

Let us briefly outline some key features of QMC needed in the present calculation. For 
a bag centred at position f in a nucleus (with the coordinate origin taken at the center of the 
nucleus), the Dirac equations for the quarks in the nucleon bag are given by |23], 24 1 



ij-d x - K - V*(f)) ~ 7° (W) ± l -V p 








The mean-field potentials at the centre of the bag are defined by V a (r) = g$a(r), V u {f) 
and V p {r) = glb(r), with g%,g% and g q p being, respectively, the corresponding quark and meson- 
field coupling constants. (Note that we have neglected the variation of the scalar and vector 
mean-fields inside the nucleon bag due to its finite size p3||.) For 3 He, which is too light for 
a simple shell model description, we use the empirical density distributions [^j| and calculate 
mean-field potentials in the nucleus using local density approximation. However, for 6 Li we 
calculate the mean-field potentials and proton and neutron density distributions in the QMC 
model self-consistently, by solving Eqs. (23) - (30) of Ref. |2T 
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The effect of the mean field potentials on the internal structure of the nucleon can be 
totally absorbed into the normalization constant, the quark eigenenergy, a small change in the 
bag radius compared with free space and the relative renormalization of the lower component 
of the Dirac spinor: 

where r' is the quark coordinate in the bag and 

E q = tt + R(y w {r)± l -V p {r)) for f^J quarks, (3) 
iV" 2 = 2R 3 j*(x)[n(n - 1) + Rm*/2}/x 2 , (4) 



P = J(n-Rm* Q )/(n + Rm*), (5) 



f2 = yx 2 + (Rm*) 2 , m* = m q — g q a a(r). (6) 

Once we have the quark wave functions in medium, we can calculate the corresponding 
parton distributions. All that are required are suitable modifications of the free space expres- 
sions, which we now summarise. 

In the case where the proton is described by the MIT bag model with just three valence 
quarks in the ls-state, the dominant contribution to the twist-2 quark distribution is given by 
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Qv(x) = j^y 3 J -x) -p+)\MPn)\ 2 . (7) 

This contribution comes from the case where the spectators to the hard collision are two valence 
quarks, again in the ls-state of the bag, which can form a scalar (total spin S=0) or a vector 
(S=l) system. The factors come from the Peierls-Yoccoz momentum projection which is 
used to build translationally invariant states. The integration is over the momentum of the 
diquark spectator system in the intermediate state, ip + is the plus component of the quark 
wave function in momentum space, and in Eq.(^) and in what follows x is the standard Bjorken 



scaling variable. More details concerning Eq.(0) can be found in Ref. |20 . 

After integration over the transverse components of the momentum of the diquark pair, 
we have: 

In this paper, we will be working with the bag model wave function for a massless quark, which 
is given by Eq.(||) with V a — — V p — for the free space case. In its simplest form, the bag 
calculation has just a few free parameters: the bag radius (for a free nucleon), the scalar (M s ) 
and vector (M v ) diquark masses, and the starting scale, fi, for the QCD evolution (at which 
the bag model is supposed to best represent the non-perturbative structure of the nucleon). 
They are fixed by fitting the valence distribution to the respective existing parametrizations for 
the experimental data. This procedure has been quite successful in the past |2(], pi]] , showing 
that the bag model is able to describe not only the unpolarized valence sector, but also has 
very good predictive power for the polarized sector p0[ . In this work, we use /x 2 = 0.1 GeV 2 , 



R = 0.8 fm, M s = 700 MeV and M v = 900 MeV, which are the values for the parameters 
which give a good fit to the MRSA parametrization of the valence distribution. 

In calculating the twist-2 parton distributions in the bound nucleon one must omit any 
interaction between the struck quark and the nuclear medium. However, the bound nucleon 
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itself and the pair of spectator valence quarks in the struck nucleon still feel the mean scalar and 



vector potentials pifl . That is, in- medium the a and lu coupling at the quark level change the 
mass and energy of the nucleon and the zero component of the momentum of the intermediate 
state in the following way: 

M -> M* + 3K? 

V l -> (p° n y + 2V«. (9) 

Note that the Bjorken variable, x, is defined in terms of the free nucleon mass, M, so that 
Mx is actually mass independent. This is an important observation as it directly effects the 
integration region. 
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Figure 1: The valence quark contributions to the proton and neutron structure functions cal- 
culated in free space using Eq. (||) and in 3 He using Eq. (|i4|). 

Incorporating these changes, the quark distribution for a nucleon at rest in the medium 
is given by: 

jn) (x) _My_f°° , \UVn)\\ j { )? { ) 

^ ' 2(M*+V%-Mx) WdK J\ 
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with 



y 



M* + 3V U 



M M 

There is a second important observation to be made regarding the in-medium calculations. This 
concerns the p meson coupling to the quarks fl23"| , |24| , ^8 



M -> M* + 3V* + \t£V* and p° n 
and quark isospin matrices and V 



N 



In this case, Eqs.fljj) are modified to 
(p°)* + 2VJ + E?=i \rl(j)V^, where r 3 are the nucleon 
Vp 1 is the rho meson potential. We note that, with the 
inclusion of the p meson, these expressions imply different integration minima for the u and 
d quark distributions in the bound proton. Moreover, in a nucleus with N ^ Z, the valence 
distributions are no longer charge symmetric between the proton and the neutron, meaning 
that the substitutions d%(x) = u™(x) and v%(x) = d%(x) may be less accurate in-medium than 
in free space. This is a model independent phenomenon, and its origin, unlike the free space 
investigations of isospin breaking in the valence distributions [53], is not in the quark masses. 



In practice, the mean field potential Vjj is quite a bit smaller than those associated with a and 
uj and we shall drop it in this first investigation. 
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Figure 2: The ratio of free to 3 He polarized, valence quark distributions in the proton for the 
u and d quarks. 

Although correct, Eq.([TT]) is not yet suitable for the in medium calculation for the fol- 
lowing reason. When calculating the quark distributions in a nucleon, we express them as a 
function of x, which is a fraction of a chosen unit of momentum (in this case the mass of the 
proton, M). If we now want to calculate the quark distribution of a nucleon which itself has a 
certain momentum distribution inside a nucleus, then this quark distribution will be a function 
of the momentum fraction y carried by the nucleon. However, our Eq. (]T0|) was derived for a 
nucleon with momentum p ^ = M* + 3K,. We therefore define a new valence quark momentum 
distribution: 

q v/fl {x) = qf\yx). (12) 

For an arbitrary nucleus, we then have the following convolution formula for the nuclear valence 
distribution at leading twist H, R3T|: 



(A) W=r4/A(y)9./)- 
Jx y 1 y 



(13) 
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where fN/A{y) is the usual momentum distribution function for the nucleons inside the nucleus] 
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This choice ensures the valence quark number (Gross Llewellyn Smith) sum rule as well 
as the appropriate momentum sum rule. The distributions presented in this paper will be 
for q v /f}(x) calculated through Eq.flL?!). We emphasise that Eq. (|T3|) is the standard nuclear 
convolution, which accounts for the kinematic corrections associated with binding and Fermi 
motion, so the ratio q v ,^{x)/q v {x) is a measure of the genuine off-shell correction to the valence 

below) can 



parton distribution. Our parametrizations for this ratio (see Eqs. (|16|) and ([II 
therefore be used in conjunction with any conventional calculation of the nucleon momentum 
distribution, f N /A(y)- 

The quark distributions in medium are not only a function of x but also a function of the 
distance, r, from the centre of the nucleus. Hence we average the quark distributions over the 
nuclear density: 

d 3 rq v m(x,r)p^(r). ( 14 ) 



X) 



Note that as we have already specified the momentum of the bound nucleon, this averaging 
procedure amounts to a semi-classical approximation. It is an interesting challenge for future 
work to improve on this point. Here Pjy(r) is the probability density for the struck nucleon 
and q v /fi(x,r) is the quark distribution calculated at the local nuclear density at r for a bag 



with the effective mass, radius and so on given by the QMC equations. Eq. (14) is valid for 
3 He, where the two protons and the neutron are in s-waves. However, 6 Li may be thought of 
as an a particle plus a neutron-proton pair, each one in a p-state. It is this pair, with deuteron 
quantum numbers, which gives the spin of the 6 Li nucleus. Thus, when calculating the 6 Li spin 
structure function we average only over the proton and neutron p-state densities (calculated 
with the help of the QMC equations). In the particular case of the proton in 6 Li, is simply 
the density of protons in the p-state. For g™ it is the neutron density in the p-state. In actual 
calculations, the structure function for the bound nucleon was calculated over a grid of values 
of r, with interval 0.04 fm, running from r = to r = 4 fm for 3 He, and from r = to r = 12 
fm for 6 Li. 
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Figure 3: The valence quark contribution to the proton spin structure function calculated in 
free space and in 6 Li. 



For the 3 He case, we use |29f pa He (r) = 0.149 exp[-0.4244 r 2 ] in Eq. (|T4|) to calculate 
the corresponding distributions. The results are shown in Fig.|l], where we show the structure 
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function, gi(x), for a free proton (neutron) in comparison with that for a proton (neutron) bound 
in 3 He. Note that we have not carried out the convolution (c.f. Eg. (|T3|) ) needed to include 
conventional binding and Fermi motion effects, so the difference between the two curves is solely 
a measure of the off-shell corrections for the bound nucleon. We see that these corrections are 
not significant in comparison with the current experimental errors. This is an important result 
when we remember that 3 He has been chosen because, from the point of view of its spin structure 
function, it is assumed to be a pure neutron target. Our result vindicates this view, since we 
get almost no correction for the neutron spin structure function. (On the other hand, our 
result says nothing new about the standard nuclear corrections associated with the kinematics 
of binding and Fermi motion, nor about the nuclear structure corrections associated with wave 
function components other than the neutron with a pair of spectator protons in the 1 So state. 
Those corrections can be estimated from standard nuclear structure calculations || j33|j.) 



0.06 




x 



Figure 4: The valence quark contribution to the neutron spin structure function calculated in 
free space and in 6 Li. Also shown is the difference g{p ee {x) - glp l (x), which is relevant for the 
extraction of g\ n {x) from the measured value of gl Lt . 

In Fig. H we show the ratios of the free to the bound polarized u and d quark, valence 
distributions. The in medium corrections are practically the same for the u and d distributions 
in the whole valence region, up to x ~ 0.45. In the large x region, the d distribution in 3 He is 
supressed relative to its free counterpart at a faster rate than the u distribution. The calculated 
ratio of the polarized, valence quark distributions for the proton in free space to that in medium 
can be parametrized in the following way: 



with the parameters: 



CLqX I CnX ( 



x 



(15) 



a u = 118.41 
b u = 18.97 
c u = 1.0758 
d u = 
e u = 0.0335 



a d = 8.964 
b d = 7.5848 
c d = 1.0515 
d d = -0.0048 
e d = 0.01086. 



(16) 
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A similar procedure is used to calculated the in-medium corrections to the proton and 
neutron spin structure function in a 6 Li target - see Figs. |3] and f|. The main differences with 
respect to the 3 He case are that the nuclear density in 6 Li is calculated self-consistently in QMC, 
as explained in the introduction, and the average is over the proton and neutron p-states only. 
The ratio of the free to bound, polarized valence quark distributions can be parametrized as: 



Aq v (x) 

with 



i q + b q x - 5 + c q x + dgX 1 - 5 + e q x 2 + f q x 2 - 5 + g q x 3 , (17) 



a u = 1.1267 


a-d = 


1.0904 


b u = -2.2012 


b d = 


-1.4308 


c u = 17.3645 


Cd = 


12.0596 


d u = -66.6193 


d d = 


-50.0592 


e u = 132.6824 


e d = 


107.9883 


f u = -131.2261 


fd = 


-116.1088 


g u = 50.7843 


9d = 


49.2399. 



(18) 

As in the case of 3 He, the off-shell corrections to be made to the deconvoluted nucleon spin 
structure functions are not large. However, in order to extract g™ (x, Q 2 ) from a 6 Li target, one 
has to subtract gi(x, Q 2 ) from the measured value of gl Ll . In this case, one can see from Fig. [| 
that the use of the free g[ (x, Q 2 ) would induce some underestimate of g™(x, Q 2 ) over the region 
0.4 < x < 0.8. The off-shell corrections to the proton spin structure functions are as much as 
6-8% of g™ itself. They will therefore become important to the determination of g™(x, Q 2 ) when 
the experimental precision gets to this level. 

Although the magnitude of the off-shell corrections calculated here is model dependent, 
we believe that the ratios should be less sensitive to the details of the model. In this sense, 
the parametrizations for Aq v (x) / Aq v ,^(x) through Eqs. flTop - (|18D are suitable for general use 
when making the nuclear deconvolution. We note also that we have verified that the ratios of 
the valence distributions presented here are approximately scale independent so that they can 
be safely used for Q 2 up to 10 GeV 2 . 

The results presented here are the first consistent calculations of off-shell corrections to 
the proton and neutron valence spin structure functions in nuclei heavier than the deuteron. 
By consistent we mean that the quark model used to calculate the parton distributions is 
the same model used to calculate the nuclear structure. All the parameters were fixed in 
order to reproduce data other than nuclear deep inelastic scattering and therefore the results 
shown here really are predictions. Most important, perhaps, is the fact that we have at our 
disposition a procedure where the off-shell corrections to the valence quark distributions can 
be systematically studied. 

Overall, we have one main prediction together with one main confirmation. First, it is 
reassuring that one can indeed identify the measured spin structure function of 3 He with the 
neutron spin structure function, without significant off-shell corrections. Second, our work indi- 
cates that future experiments using 6 Li as a target, will need to take into account the in- medium 
corrections, of the type calculated here, when extracting g"(x). 
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